A New Method for Obtaining Transient Solutions of 
Electrical Networks 

By W. P. MASON 

Synopsis: A new method for obtaining transient solutions of electrical 
networks is developed in this paper which depends upon the fact that a 
distortionless line can be made to approach as a limit all three of the circuit 
elements, resistance, inductance and capacity. The process of solution 
consists in solving for the current in a distortionless line — which is ordinarily 
a simple process — and then proceeding to the limiting case of the distortion- 
less line which approaches the element or elements of interest. Some ex- 
amples are worked out and a derivation of the Laplacian integral solution 
is given. It is interesting to note that this method gives a formal solution 
of the Laplacian integral equation. 

THE following paper sets forth a new method for obtaining the 
transient solutions of electrical networks, which it is believed 
has some advantages over other methods of solution, in that the 
operations required for solution are quite simple, and also because this 
method presents a more definite physical picture of the processes 
involved. By means of this method, the current at any time t can 
be obtained, due to an applied voltage which is zero when / is less than 
zero, and is E cos (cot + 0) when t is greater than zero. This type 
of voltage includes as a special case the applied voltage, which is 
zero when / is less than zero, and is unity when / is greater than zero, 
and hence the solutions obtained by this method reduce to the cases 
discussed by Heaviside, 1 when co and are taken equal to zero. 

This method gives directly the more compact Laplacian integral 
equation solution, first obtained by Carson, and in addition gives a 
method for solving this integral equation, if its solution is not already 
known. 

I. Method of Solutiox 

All practical schemes for solving the transient type of circuit 
problem, including the Laplacian integral equation, and the generalized 
Fourier integral solution, are made to depend on the known and easily 
determined steady state solution. This implies that all circuits which 
have the same steady state solution, have also the same transient or 
time solution. The method described in the present paper rests on 
the same basis. 

The method of solution used here depends upon the fact that the 

1 Heaviside, "Electromagnetic Theory," Volume II. 
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distortionless line can be made to approach as a limit, all three of the 
electrical elements, resistance, inductance, and capacity, and that the 
complete solution for the current in a distortionless line can be 
obtained by adding the incident current and the sum of the reflected 
currents which can occur up to the time of interest. That is the dis- 
tortionless line has a true velocity of propagation, and hence the 
current at any time will be the initial current and the sum of the 
reflections which can occur up to the time of interest. All of the 
three electrical elements, resistance, inductance, and capacity, can be 
considered as limiting cases of the distortionless line. Hence the 
process of solution consists in solving for the current in the distortion- 
less line, and then proceeding to the limiting value of the line which 
coincides with the element of interest. 



¥ ¥ 



Fig. 1-A. 
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Fig. 1-B. 
Diagrammatic representations of lines. 

A . The Distortionless Line 

Since the distortionless line is the tool by means of which problems 
are solved by this method, a brief discussion of lines 2 is given here. 
If a voltage is suddenly applied to a transmission line, the current 
at any point in the line is zero for a certain time and then begins to 
build up to its final or steady state value. If there is no distributed 
inductance in the line, the current begins to build up immediately. 

For a distortionless line, however, the current is zero for the time 
required to propagate the wave to the point of interest and then 
instantaneously reaches its steady state value. To show this let us 
consider the equations for a transmission line. A line has distributed 

2 For a more complete discussion of lines see "Transmission Circuits for Telephone 
Communication," K. S. Johnson, page 144. 
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series resistance R u and inductance L,„ and distributed shunt capacity 
C u and leakage conductance G u , as shown on Fig. 1-A, where the letter 
11 indicates the values per unit length. If i denotes the current and 
v the voltage at a distance x from one end of the line, the well known 
differential equations are 






(i) 



If we eliminate v from these two equations we have 

L u Cu ^ + (RuC u + G U L U ) f t + R U GJ =~ (2) 

Similarly, if i is eliminated, the resulting equation is the same as (2) 
with i replaced by v. Since we are dealing with simple harmonic 
forces, the current i varies as cos (co/ + 0) where a> is 2w times the 
frequency of vibration and is an arbitrary phase angle. It is 
usually more convenient to consider i as varying according to the 
time factor 

i = yw+t = |[cob (at + 9) + j sin (w/ + «)], 

where I is the maximum amplitude of i. The solution obtained on 
this assumption is called the symbolic solution, and the real solution 
is obtained from the symbolic solution by taking the real part. 
Substituting the symbolic form of i in equation (2), this equation 
reduces to 

[0)--L u C, + MR U C U + G U L U ] + R u GuJ = ~ (3) 

The solution for a line can be specified in terms of two parameters, 
the characteristic impedance and the propagation constant of the 
line. To show this we note that the solution of (3) is 

i = Ae~ rx + Be**, (4) 

where T 2 = [i? u + juL u ~][G u -f- jwCu] and A and B are constants. 
From the last of equations (1) we have 

A • 

» = ~ r m- r = r I- r ^ Ae ~ Tx ' Be ^ ^ 
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When x = 0, i = i , and v = v . From (4) and (5) solving for A and 
B we have 

vj _ fo_i_ ito/2 . R = i? Po/2 



' G u + jo) C u \ G„ + jw C u 

Substituting these values in (4) and (5), we have the equations 

, _, v sinh Tx 
i = to cosh I x , = > 

i ^, • /./?„ + jcoL H . . r, 
v = v cosh Tx — to \ \ r ■ ■ .^ sinh 1 x. 

In this equation Tx = P, the propagation constant of the line, and 



4 



l R u + juLu = z 
Gu "f" iwC„ 



the characteristic impedance of the line. If we are interested in a 
given length of line /, the parameters are 



\R + juL 



P = Tl{R + j a L){G + j»C); Z Q - VGTfcc' (?) 

where R, L, G, and C are the total distributed constants for the length 
of line considered. The characteristic impedance is the impedance 
looking into a line of infinite length as can readily be seen from either 
of equations (6) by letting x, the length, approach infinity. In this 
case cosh Tx = cosh P approaches sinh P and both approach infinity. 
Then from (6) 






/COSh P — ilio\ ry ~ 1 

. . „ Zo -> Z when x -> » 

\ sinh P J 



since i/i can never be larger than 1. 

The physical significance of the propagation constant is that e~ p 
represents the ratio of currents or voltages at the two ends of the line 
when the line is connected to an infinite line of the same characteristic 
impedance. To show this, suppose we terminate the section con- 
sidered in an infinite line, which as we have seen above has an im- 
pedance Zo- Then in equations (6), v = »i the output voltage and 



TRANSIENT SOLUTIONS OF ELECTRICAL NETWORKS 113 

i = ii the output current. We let Vi = iiZ . Eliminating either v 
or to, we have 

V\ = v e~ p or i\ = ioe~ p . (8) 

In the following work it is necessary to know the impedance of a 
short circuited line and that of an open circuited line. For the short 
circuited line, the voltage at the far end is zero, so putting v = in 
the last of equation (6), we have 

Vo/io = Z tanh P. (9) 

For the open circuited line we put i = 0, obtaining 

Vo/io = Z coth P. (10) 

So far we have discussed the general transmission line. For the 
distortionless line there exists the relation 

r-i- <»> 

Substituting this relation in equation (7), these parameters reduce to 

Zo=-Ro = J! = vl and P = JRG+Ju^LC = A+juD. (12) 

This equation shows that the characteristic impedance becomes a 
resistance Rq, while the propagation constant equals a real constant A 
plus ju times the constant D. To show how wave transmission 
takes place in an infinite line, these values are substituted in equation 
(8), giving 

To find the real solution, we take the real part of this symbolic solution, 

obtaining 

Vl = v a e^ l+B) e- u+ ' aD) = v Q e- A e'l^-D)+e ]i 

or, taking the real part, 

Vi = v Q e- A cos [«(/ - D) + 0] 

and (13) 

/, = i e~ A cos [w(/ - D) + 0], 

where the clash over z>o and «o indicate the maximum amplitude of 
these quantities. Since v = £0 cos (cot + 0), these equations show 
that either V\ or /, has the same form as v or i respectively, attenuated 
by a factor e~ A and delayed in time by an amount D. 
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B. Condition for Obtaining Lumped Constants from a Distortionless Line 
In a distortionless line there is only one necessary relation between 
the constants, equation (11). Hence, we are at liberty to vary the 
constants subject only to this relation. In the following work we wish 
to make the distortionless line degenerate into resistances, inductances, 
and capacities, or combinations of these. 

For example, suppose that we wish to obtain a resistance from a 
distortionless line. To obtain this we take a short circuited line as 
shown on Fig. 1-A and let R be finite, L -> 0, G -^ 0, and C -> 2 in 
order to satisfy equation (11). The shunt elements will all vanish 
and the series inductance disappears, leaving only the series resistance 
in the line. Since the line is short circuited, the line degenerates into 
a resistance. The line parameters for this case become 



R, = A IE = Jk -> oo 1/2 ; P = jRG + j<* VLC -> (0" 2 + >0»«) 
Y G Y C 



and 



Ro P = Jp (,[RG) + J? {jujLC) -> R. 



(14) 



There are three combinations of lumped elements which can be 
obtained from the short circuited line and three combinations which 
can be obtained from the open circuited line. These are listed in the 
following table, together with the resulting line parameters. 

Limiting Cases of the Distortionless Line 
A . Limiting Values with Short Circuited Line 



Equa- 


Assumed Line 
Constants 


Resulting Line Parameters 


Resulting Lumped Element 




R 


L 


G 


C 


Ro 


P 


RaP 




(14) 


Finite 








2 


oo 1/2 


Ql/2 + j w 3 / 2 


R 


Resistance • — WW"»-« 


(15) 





Finite 


2 





00 1ft 


Q3/2 _|_ jaQllS 


joiL 


Inductance •^TRnnT^ 


(16) 


Finite 


Finite 








00 1/2 


Ql/2 _|_ j w Q 1/2 


R+joiL 


Resistance y\y\^_TK^r\^ 

and 
Inductance 



The first three cases result from the suppression of the shunt 
elements in the short circuited line, and the line parameters are 
characterized by R -> » ; P -> 0. The last three cases result from 
the suppression of the series elements in the open circuited line, and 
are characterized by Ro — > ; P — > 0. 



TRANSIENT SOLUTIONS OF ELECTRICAL NETWORKS 115 
B. Limiting Values with Open Circuited Line 



Equa- 


Assumed Line 
Constants 


Resulting Line Parameters 


Resulting Lumped Element 




R 


L 


G 


C 


Ro 


P 


RolP 




(17) 





2 


Finite 





i/2 


Ql/2 + j u 3 ' 2 


l/G 


Resis- 
tance *wvw-» 


(18) 


2 








Finite 


01/2 


Q3/2 + juQUl 


IfjuC 


n . ,| 




(19) 








Finite 


Finite 


Ql/2 


Ql/2 + jtaQtH 


l/(G-h7«0 


Resis- #_{ 1 1 (__. 
tanceand ^VW^ 
Capacity 



C. Solution for a Resistance and Inductance in Series 

As a first example of a transient solution obtained by this method 
let us consider the case of an inductance and resistance in series with 
a source of alternating voltage. To solve this problem, consider the 
case of a voltage in series with a distortionless line, short circuited, 
as shown on Fig. 1-A. The current immediately flowing on applica- 
tion of the voltage will be i where 



to = 



E 

Ra 



(20) 



This current is transmitted down the line and completely reflected 
at the far end, returning to the near end. The first reflected current 
entering the generator is then 



t\ = t e~ 



(21) 



where P is the propagation constant of the line. Upon reaching the 
near end, the current is completely reflected in the same phase and 
again enters the line. At the end of the first reflection, the current 
entering the line is 

i = * (1 + 2e~ 2p ). (22) 

After (n — 1) reflections and passages through the line, the current is 

i = j [l + 2e~ 2P + 2e~ ip + • • • + 2e-^ n ~ 1 ^ 

(23) 



-«KV^)-'] 



Now the time at which the nth reflection occurs will be 

/ = n(2D), 
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where D is the time of delay in passing the network once. For a dis- 
tortionless line 

D = «$LC. (24) 

Hence, we can replace n by 

-arsra- (25) 

Since D -» and n — > °° , the time scale becomes continuous. In 
equation (23) we insert the values given in (16) and (25), appropriate 
to the limiting case considered here, namely 

n R + jtaL D t tR 

and note that 2P -> so that e~ 2P -> 1 - 2P\ then 



^ r 2(1 — e -tRoiuui+juL)iR }') j 

1 -*R [ 1 - 1 + 2(R+jcoL)IR ~ 1 \ 



(26) 



r 1 - e- tiRiL+ ^ _ J_ 1 
~ L ^ + >L i?o J ' 

But i?o -* °° and hence the solution is 

[J _ e -t(RlL+j") "1 
R + j«L J' 

This is the symbolic or complex algebra solution of the equation 

L%+Ri = E. (27) 

at 

In general it is desirable to obtain the current due to an applied 
voltage of the form 

E = E cos («/ + 9). 

This solution can be obtained directly from the symbolic solution 
given in (26) by taking the real part. The result is 



* = E, 



cos (<at + - g) ~ cos (fl - <p)e- tRlL "j , 2g s 

V-R 2 + co 2 L 2 J 



It will be found that (28) is a solution of (27) for an applied voltage 
£ cos (orf+ 0). 
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D. General Method for Determining Reflections 
The method for obtaining the successive current reflections of the 
line given in the preceding section, is laborious to carry out in com- 
plicated cases and hence it is desirable to obtain a simple method for 
determining the reflections. Such a method is the expansion of the 
expression for the impedance by the binomial theorem in order to 
get the successive reflections. In the above example the current i is 

* = E/(R+ja>L). 

We note that the expression for the impedance is approached by that 
of a short circuited line when the R and L of the line are finite and the 
capacity and leakance approach zero. Hence 

E E _ £(! + e~- p ) 

R + ja>L ~* R tanh P ~ R (l - e~ 2p ) 



1 



Now the expansion of 
1 



1 - e~- p 
Hence 



= 1 + e~ 2p + e~ ip + 



i -»-f (1 + e-* p )(l + e~"- p + e~* p + • • •) 



R 



= ■§■ [1 + 2e~ ip + 2e~ 4p + • • • + 2e-^»-» p +•■•], 

which is the expression for the reflections given by equation (23). 

In all the following problems it will be found that a similar process 
for obtaining the reflections can be followed. It is evident that any 
method which gives an expression for the current in the form 

i = £[a + air*** + <hfi-*° D + h a n e- n -'^ D + ■ ■ • ] (29) 

will give the reflections, for if we take the real part of this expression 
we have 

* = £ [>o cos (ut + 6) + fli cos [«(/ - 2D) + 0] + • • • 

+ a n cos [«(/ - 2nD) + 0] + • • •]• 

Each term represents a current which adds to the original current 
after a time of delay 2D, AD, ■■ • • 2nD, and hence the wth terms repre- 
sents the wth reflection. Therefore any method, such as the above, 
which gives the current in the form of equation (29), will give the 
reflections. 
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E. Simpler Form for Replacing an Impedance 
In the preceding section, the transient solution of an inductance 
and resistance in series was obtained by replacing the impedance 
R + jaL by the expression 

R tanh P, where R P = R + jwL and R -> « ; P -> 0. 

Tanh P has a numerator and a denominator both of which must be 
expanded in order to obtain the reflections. If a single term can be 
used, the expansion becomes simpler. In order to effect such a 
simplification, it is necessary to find a physical structure, which has 
only one term in its impedance expression and which approaches a 
resistance and inductance as a limit. 



Ro' 




R0- 00 p-o 

Ro P = R+jOJL 



Fig. 2 — Short circuited line and shunt resistance. 

Such a structure is shown on Fig. 2. It consists of a short circuited 
line shunted by a resistance R . The current into the combination is 



i = 



E 



E 



Ro X R tanh P P„(l - or**) 



(30) 



R + Ro tanh P 2 

If now in the short circuited line, we let R and L be finite and G — > 0, 
C — > 0, the combination obviously reduces to a resistance and induc- 
tance in series, since the infinite shunt will not affect the result. 
Hence the replacement of a resistance and inductance in the equation 

• E 

1 R+jo>L 

by the expression in (30), is justified. 
The solution of (30) is gotten by expanding the expression and is 






-2(n-l)P 



+•••]= 



2£[1 - e" 2nP ] 



Roll - e-™~\ 

Upon substituting in the values PoP = R + jo>L, n = t/2D, and 
letting R —> °° ; P — > 0, we have 

_ e -t[RIL+ia] ~] 



i-s[l 



R + juL 



in agreement with equation (26). 
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Similarly, when we have the expression 

1 



we can replace it by 
2R Q 



1 - e~ 2p 



, Ro 
where -=r = 



G + jcoC 

_1 

P ~ G + jaC 



and i? ->0, P -» 0. (31) 



The structure which gives the impedance in (31) is an open circuited 
line in series with a resistance R . The impedance of the combination 



is 



R + Ro coth P = R [ 1 + j4"pS] = i"^ ' 

If then the series impedances of the line approach zero, R — > and 
the impedance of the combination approaches 

_1 

G + jaC 

F. Solution for a Resistance and Capacity in Series 

As a second example let us obtain the solution of a resistance and 
capacity in series. To obtain the solution we solve the case of a 



R — O P—o 

R °/p^yjLoc 



Fig. 3 — Open circuited line and series resistance. 

resistance in series with an open circuited line as shown on Fig. 3. 
The steady state solution for the current in this circuit is 




E 



Replacing -i—r, by 
JcoC 

there results 
E 



9 7? 

^ , and substituting in the above equation 

1 — e ^ 



R + 



2R G 



when Ro -> and P -> 0, and -^ -^>— 



P juC 



1 - e~ 2P 
in accordance with equation (18). 
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After some rearrangements this can be put into the form 

E [ R e-"(l + e-*n I 
* " i? + R I * i? 1 - e- 2 <™-> J ' ^ ; 



where 

g-2F _ 



Expanding equation (32) in the form of a series, there results 

i - E |~ -R+3i?o _ Rj> / 2R+2R \ , 2(F+P ), -4(f+p) . \ 1 

Summing up n terms of this series, we have 

E r R + 3R R 1 2R + 2.Ro \ / 1 - e- 2 "'™ 1 ) 



« = 



i? + i?o -K + 2R R\R+2R /\1 - e-^+ p > 



(S3) 



Since in the above expression Rn — > 0, we can obtain the value of F 
by writing the first terms of the expansion for the exponential 



Hence 



r* = 1 - 2 F + S3 + • • • = * = 1 - ^ + 



-4° 



If now in equation (33) we proceed to the limit, letting 

i?o->0; 
there results the equation 



S ->0 : P-»0; ff-j^ji ,-jJ, 






1 _ e -t(URC+M 1 

'- i-Uot I" (34) 



This equation is the symbolic solution of the integral equation 



Ri + i fidt = E e^'+ 6 K 



If we wish the solution corresponding to the impressed voltage, 
Eo cos (cot + 6), we take the real part of (34), obtaining 



_, cos (<at + 6 + 8) - [sin (0 - 5) tan «]«-" 

= En . 



, .,,, .... ,,„ , ... ,,,:. „ ,f- llRC 

i = Eo , =- 

Vie 2 + i/co 2 c 2 

where 

tan 5 ■ ^k ■ 
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II. Solution for M Sections of All- Pass Lattice Network 

The process for solving any type of problem is to replace any 
resistance and inductance in series by l?o,(l — £ -2/>I )/2. and any 
capacity and leakance in parallel by (1 — e~- p -)/2Ro 2 , where Ro t — > °° ; 
R , -> 0; R Ql Pi = R + jwL; R^IP* = l/(G+j(oC). ' 

The next problem considered here is the solution for any number of 
sections of all-pass lattice type network s as shown on Figure 4. 



*=V^c 




r=Vl^ 



Fig. 4 — Sections of all pass lattice network. 



These networks have the property of passing all frequencies without 
attenuation, and they are much used as phase equalizers. 

The steady state equation for the current at the end of the first 
section, when this section is terminated at each end by a resistance 
R = VZ1Z2, is 

e r Vz;- Vzx ] 



h = 



The current flowing out of the rath section of such a structure takes 
the form 

E rVzl- VFi' 



. e r vzo - jgj T 

*" 2i?|_Vz 2 +VzJ 



(35) 



In the structure considered Z\ = juL; Z 2 = ljj**C, and *lLfC - R- 
In accordance with section I-B, we replace the inductance by a short 
circuited line, and the capacity by an open circuited line. For the 
first line, in the limiting case, we have by equation (15), 

R -> =c ; RoP = juL. 

For the second line, by equation (18), we have 

P jut 

There is no loss of generality if we take the propagation constants for 
3 See for example B. S. T. J., July 1928, page 510. 
9 
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the two lines equal so that 



R 



L 



^ X RoP = j<*L X -™ - £ = R 



P 



jcoC C 



Hence R = R 2 /'Ro- Substituting these values in equation (35), we 
have 



2R[l 



-[i?o/2i?](l -e-"-nY 



I 



(36) 



+ [i?o/2^](l - e-* p ) 
After some simple rearrangements, equation (36) takes the form 

im = 2R[R~ \l-e-™<«^)~ l \ ' (3?) 

If m equals 1, the solution can be obtained exactly as discussed in the 
first example in section (1), and it is 



E f 1 - jwjLC - 2e- tl ^ LC) + ]u] 



H 2R\_ 1+>a/LC 



] 



The solution for m sections of lattice network is discussed in the 
Appendix, and it is there shown that the solution can be written in 
the form 



E_ 
2R 



L (m-l)! 1+joWLC 



\l+ja>^LC/ 

(-±=y- 2 

■ \4LC) ( 2- n 

"•" (m - 2) ! \ (i + jaJLCy 1 + jwVL C / 

j r q-jWIC)- _ n / , 1 2L_Y| 

"*" ' vIcL(i+>VLC)- 1 \ l+i«VLC/J 



wj*- 1 \ 



+ 



Vl+i 



-jodLCy _(_!), 



joi-yLCj 



(38) 



Equation (38) represents the symbolic or complex algebra solution 
for the current at the end of the wth section of a lattice network as 
shown on Fig. 4. It is usually desirable to obtain the current due 
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to an applied voltage E cos (ut + 0). This can be obtained from 
equation (38) by taking the real part of the equation and rejecting 
the imaginary part. The process of doing this is simple and the 
result obtained is 



t m = 



2R 



+ 



cos (ut+6 - 2m<p) - 2e-«'< LC > 

(vzlO 



/ _2/_ \ m ~ l 



cos (0 — cp) COS tp 



(m - 2) ! 
/ JU_ \"'- 3 

Wlc 1 



(m- 1)! 

£2 cos 2 <p cos (0 — 2^) — m cos <pcos (0 — ip)~\ 



*" C J, 4cos 3 <pcos (0 - 3<p) - 2mcos 2 pcos (0 - 2?) 
(m - 3) ! L 



. mint — 1) /fl x I i 

H ^-yj -cos <p cos (0 - <p) + 



(39) 



where tan ip = uyLC. 

For example, the solutions for one and two section networks take 
the form 

fi = ^[cos (a./ + - 2<p) - 2erWfi& cos ^ C os (0 - p)] 

* " ik \ C0S (co/ + d " 4<p) ~ 2e " (l/VZ5) [ 2 cos2 * cos (* ~ 2 *) ( 4 °) 



*2 



2/ 
— 2 cos <? cos (0 — tp) + 77f7; cos ^ — W cos * 



]]■ 



It appeared desirable to obtain some numerical calculations on the 
building up of current in this type of network. This calculation has 
been carried through for two section, four section, and six section 
networks. The current has been calculated for a constant voltage, 
for an alternating voltage whose frequency is the resonant frequency 
of the network, and for one whose frequency is twice the resonant 
frequency of the network. The current building up for a constant 
applied voltage is shown for the three networks on Fig. 5. The 
current building up in a two section network, in which an alternating 
voltage whose frequency equals the resonant frequency of the net- 
work, is shown on Fig. 6. The steady state and the transient terms 
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are shown separately in the dotted lines, and the complete solution in 
the full line. The applied voltage is of the form E cos ut and hence 
6 is taken as zero in equation (39). Similarly, curves for two, four, 
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VALUE OF TIME X CRITICAL FREQUENCY = tfc 

Fig. 5 — Current resulting from the application of a constant voltage, E, on several 
sections of lattice network. The current plotted is the current in the termination of 
the network. 
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Fig. 6 — Current resulting from the application of an alternating voltage, 
E = E cos w c t, on a two section lattice network. The current plotted is the current 
in the termination of the network. The frequency of the applied voltage is the 
resonant frequency, f c , of the network. 

and six sections are shown on Fig. 7 and Fig. 8. Fig. 7 shows the 
transient terms and Fig. 8 the complete solution. Fig. 9 and Fig. 10 
show similar curves for a frequency twice the resonant frequency of 
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the network. In addition, the solution for infinite frequency is 
readily obtained from (39) since for this frequency <p = 90°. Then 

E 

i = ^-^ Tcos (cot 4" 9 — wtt)"]. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 OS 09 1.0 1.1 1-2 1-3 '•* 15 
VALUE OF TIME X CRITICAL FREQUENCY ■ tf c 



1.6 1.9 



Fig. 7 — Transient current resulting trom the application of an alternating voltage, 
E = E cos to c t, for several sections of lattice network. The current plotted is the 
current in the termination of the network. The frequency of the applied voltage is 
the resonant frequency, f c , of the network. 



III. Laplacian Infinite Integral Equation and Its Formal 

Solution 

The solution of circuit problems by means of the Laplacian integral 
equation has been used by Carson * to a large extent. It is interesting 
to note that this integral form can be derived in a simple manner by 
means of this expansion method, and that this method provides a 
means for solving the Laplacian integral equation. 

Any impedance Z is made up of resistances, inductances and capaci- 
ties, and hence the current i can be represented by a series 



i - 1 = £[>„ + air*® + a*-** + 



■ffl n «-' 2 -n •••]• (4i) 



The interpretation of this expansion from a physical standpoint 
is that the current is Ea , for the first interval of time 2D, 
E[_ao + aie~ iimD '] for the next interval of time 2D, etc. Hence at 
the time / = n(2D), the current i will be given by the sum of n terms 

4 See "Electric Circuit Theory and the Operational Calculus," B. S. T. J., 
October 1925, and following. 
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of this series. We can, therefore, express the current i at the time 
/ by the integral 

• = E T C a(t)er*"dt + a 1 , (42) 

where the value of ait) for any interval of time (« — l)2D to n(2D) 
is the constant of the above series a„_i divided by 2D. The value 
of this integral for an infinite time must reduce to the steady state 
value of i = E/Z, hence 
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[f- 



{t)e- iut dt + a 
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Fig. 8 — Current resulting from the application of an alternating voltage, 
E = E cos u c t, on several sections of lattice network. The current plotted is the 
current in the termination of the network. The frequency of the applied voltage 
is the resonant frequency, f c , of the network. 

Cancelling out the common factor E, we have the infinite integral 
equation 

z fe-[jf-cor*«+*]- («) 

The physical interpretation of the quantity a{t) is readily obtained 
by reference to equation (42). If we set co = and E = 1 in this 
equation we have 

i = f a[t)dt + fl = f a(j)dt + ^(0), 
Jo Jo 
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where h(0) is a constant denoting the current when / is zero. Now i 
at any time t is the indicial or direct current admittance, designated 
by h(t), hence a{t) is 

<t) -|c*(0). 



The infinite integral equation (43) takes the form 



1 



ZO) 



tr 



| (h(t))e->"'dt + A(0) 



(44) 



This integral equation does not have quite the same form as Carson's 
integral equation but can be readily put into that form by means of 
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VALUE OF TIME X CRITICAL FREQUENCY = +f c 

Fig. 9 — Transient current resulting from the application of an alternating voltage, 
E = Eo cos 2a} c t, on several sections of lattice network. The current plotted is the 
current in the termination of the network. The frequency of the applied voltage is 
twice the resonant frequency, f c , of the network. 

Borel's theorem 5 which is given below. Suppose that 1/Z' and 1/Z" 
are two admittances, which when multiplied together give the admit- 
tance l/H. The admittances 1/Z' and 1/Z" have the expansions 



-^ = [ fl0 " -f ai "e-°-' uD + a 2 "e-^ aD + 
6 See B. S. T. J., October 1925, page 722. 
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The expansion \jH, then takes the form 

Ti = Wz Ff = Lao + aie ~ 23UD + • " * + a » e ~ MD + ' ' "3 

= [>oV + [aoV + aWy**™ H 

+ [>oW + • • ■ + a fc 'a„-*" + • • ■ + oWy**" + ■•■]■ 

We have then the relation 

a n = [flj'fln" + • • • + a k 'a n -k" + • • • + dnWl- 

Now a n " is the value of a" when £ = «(2D), hence the above relation 
can be put into the form of an integral 

a{t) = | a"{r)a'(t - r)dr, 
Jo 

where r = Jfe/2Z>, and the complete relation is 

= fT fa"(r)a'(t - r)dA e-> at dt + a W. 

Suppose now that we let Z' = jw; Z" = Z{j<x>). We know from 
Heaviside's rule and from Section I that 1/ju has the direct current 
or indicial admittance solution 

i- = t = h(t). 
Hence 

a '(t) = j t h'{t) =j t t=l and a ° = °» 

and the infinite integral equation 

ti =iJw) - r [X'°" (T)a ' (/ - T)d *] ^ + aM ' 

takes the form 
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Dropping the primes, we have the Laplacian integral equation 



1 



jo)Z(jco) 



/•BO 

Jo 



h(t)e-> ut dt. 



(45) 



Hence (43) is equivalent to Carson's integral equation, if (ju>) is 
replaced by p. 

It will be noted that in deriving this equation use is made only of 
the general form of the expansion of admittances. For particular 
admittances, the values of the a's in equation (43) or the h's in equa- 
tion (45) can be derived directly from an expansion of the admittance 
function, as shown in the foregoing work. Hence, if the solution of 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 08 0.9 1.0 II 1.2 1.3 1.4 15 1.6 1.7 1.8 1.9 
VALUE OF TIME X CRITICAL FREQUENCY = tfc 

Fig. 10 — Current resulting from the application of an alternating voltage, 
E = E cos 2u c l, on several sections of lattice network. The current plotted is the 
current in the termination of the network. The frequency of the applied voltage is 
twice the resonant frequency, f c , of the network. 

the integral equation is not known from a table of integrals, one method 
for obtaining its solution is the expansion method developed above. 
This method may then have some application as a method for solving 
integral equations. 

A . Illustrative Example 
As an illustration of the use of this method in solving integral 
equations we will consider the equation 



1 == r a(t)e->»'dt. 

V(»- + 2X/'w Jo 



(46) 
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The expression on the left can be written 



1 (47) 



A^coV2X -f~ j& 

Noting that the square of the first factor has the form of an inductance 
and the second the form of a resistance and inductance in series we 
replace 

1 — e~ 2Pl \ 

where R , -> *> ; Pi -> 



jta -> Ro, I 

Ro.Pi = > 

/ 1 — e~ 2Pl \ 
and 2X + j<a -> R 0l I ^ ) where R a* ~* °° ! p * "* °» and 

i?o 2 -P2 = 2X + jo. We note that Pi has the form juD where R 0l D = 1 , 
while P 2 has the form A + joiD where i?o 2 ^4 = 2X and R ,D = 1. 
Substituting these values in (47), we have 

1 



Expanding these two factors by the binomial theorem, we have 



2 



VP^L 1 + 2 g + 2! 

(2n) !e-« 2nw ^ 



+ 



] 



■+i 



2 2 "(w!) 2 

(2m) !e-*»u-w«M 



X 1 + s e- 2 "+*"» + • ■ • + 



2 2 »(w!) 2 



At this point we make use of Stirling's theorems on factorials which 
states that when K is large 

K\ K 



K\ = Ijj 42ttK. 
The typical terms of the above expression become 

(2«) Ig-^q-w _ V g / ^ e-'< i,wD > 

22 "( w! ) 2 2-^V n 2.« ^ « 
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Inserting this value in the above expression and multiplying, there 
results 



( 



1 + e-'^^ll/2e- 2A + 1/2] + h e _,(2 '" 



'/)) 



-,-2(n-l)4 p -2(n-k)A 

+ ■■■ + .,! ■.,.,-■■+••• W 



■<[jm 2V7r(w — 1) V7r(w — &) V^ 



■a—] 



■VttW, 

Since the value of a{t) is given by the factor multiplying the term 
e -i2no>D divided by 2D we can write 

2 Ve~ inA . e -2(n-DA , 

a(t) = , — — = 

'™ 0l R 0l 2D LVttj 



-jRofio^D |_ Vttw 2Vtt(w - 1) 

7TWJ 



»-2(n-fc)A 

•\7r(w — k)-\irk ^l^ 

This expression can be written as the sum 

i=» /,— 2(n— k)A 

t=oH(n — k)k 

We introduce now the value w = //2D and define a new variable t 
by & = t/2Z). Inserting these values in the above summation and 
noting that A/D = 2\/R Jl/R 0l = 2X, we have 

fc=J. e U-r)iK 2 D 

But 2£> = dr, the element of time, so that the summation can be 
written as the integral 

The value of the integral is Te Xl I (\t), where Io(kt) is the Bessel's 
functions Jo(i\t)> when i = -V — 1. To show this we can expand 
the exponential and integrate the series giving 

'[*+»+W +i W + "]- 

This can be recognized as the series expansion of e Xl I (\t). Hence 
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the value of (51) is 

«(/) = r^JoCXfl, (52) 

which is the solution of the integral equation. 

IV. Other Types of Boundary Conditions 

The solutions obtained before are all on the assumption that no 
energy exists in the network before the voltage is applied. Other 
types of boundary conditions are sometimes desirable, but these can 
all be derived from the above solutions. 

The next most important case is the case where the network has 
come to its equilibrium value and the voltage is suddenly taken off. 
This condition is the same as would result if a negative voltage E 
were suddenly applied to the circuit, and hence the solution is the 
steady state value of the current minus the current which flows on 
application of the voltage E. 

Another type of boundary condition sometimes occurring is the one 
where energy exists in the network when / is zero. This may be 
taken iiccount of by assuming that the voltage is applied before t 
equals zero. To take account of this condition analytically, examine 
the expansion 

i = | - £[a + a ie -^ D + a*-**"* + h a n e~°-"> aD + • • •]• 

According to the above assumption, the voltage is applied when 
/ = — to, hence for n we substitute 

_ J_, _/o_ 
n 2D^~2D' 

The above series is then replaced by the integral 

i = eC a(t+ to)e-'^ l+ 'o ) dt. (53) 

J -to 

Another boundary condition of interest occurs when the voltage 
is taken off before an equilibrium value has been reached. If we 
count time as starting when the voltage is taken off, or what amounts 
to the same thing, when a negative voltage is applied, the symbolic 
solution takes the form 

i = El f a(t + to)e-'^ l+, o ) dt - f a(t)e->*"dQ. (54) 

J -to Jo 
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Appendix 
1 



The expression 



m 



4R 

R V 1 - g-2 (*/*>+« 



- 1 



can be expanded into the form 

T 1* / 1 \ 1 '" _ [ 4R I 1 



(,Rim+p) 



+ 



+ (- I)'"- (55) 



and hence the general solution depends only on the solution of the 
general form 

47?/ 1 



m 



R Q V 1 - e-2<*'A>+« 



(56) 



If equation (56) is expanded by the binomial theorem, there results 
the expression 



**y [ i + ^w. + *&+!>,-<««* + 



+ 



(m + « - l)!e-2«c«/«o+p> 



w!(w — 1) ! 

For any value of m, we can write the typical term of (57) as 

( m + n - l\»'+»- 1 
[m + n - 1) ! 
nl(m - 1)! 



(57) 



V2tt(w + » - 1) 



I2irn(m — 1)! 



( m+ ;;- 1 )"( W + >/-i)»-^ 



(58) 
(m — 1) !V« 

by Stirling's theorem on factorials. Now n for any finite value of time 
approaches infinity, while m for any finite term in the series is finite. 
Hence (58) can be written as 



e -(w-l) ( i _|_ 



m — 1\" 



(w - 1) ! 

m 1 \ " 

The limit of ( 1 -\ I as n -> « is e*'"- 1 ). 6 Hence 



(m + n - 1)! 



«!(w - 1)! (ra - 1)! 

6 See "Probability and Its Engineering Uses," T. C. Fry, page 107. 
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The value of (57), then reduces to 

i(RlRo+P) 






(59) 



If now we substitute n = t/2D, P = jcoD; R^P = juL (59) reduces 
to the integral 



4^ \ m J nt tm-i e -iai VZc)+/«] « 



(OT - 1) ! 

2_ 
"LC, 



X t — — tt7 r t^e-wi^+Mdt. (60) 
(« - 1)1 Jo 

If we integrate (60) by parts, successively, there results the series 



2 m . 



1 



1 +jw4LC 

+ 



_ />-'(!/ VZc+i") 



\VZZ7 



(1+ JutILC)™ 

f \ m—1 



(1+jwVZC)- 1 



+ 



Vlc 



(w- l)!(l+icoVLC). 



(61) 



The complete solution of (55) is then obtained by adding terms of 
the kind given in (61). For example the steady state term is given 
by the series 



w2 m -' 



mijn-V) 



2! 



=— - ■+■■• + (-!)■ 



(1 + >VLC) m (1 + jwVrC)"- 1 (1 + co^LC) m ~ 2 

2 m - w2— »(1 + jcWIC) + m( - m 2 ~ !J 2"'- 2 (l + ju-yJLC)* 

m +_••• + (- !)-(!+ jwVEC)' 

(1 + jco^LC) m 

This is readily seen to be the binomial expansion of 



f 2 - (l+jt,rfLC) ~\ m _ f l -jctfVLC I" 
L 1+jaWZC J Ll+jcoVZcJ 



(62) 



The other terms given in equation (38) follow directly by addition 
and reference to equations (55) and (61). 



